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We use Fermi liquid theory to study the mechanical impedance of '^He-*He mixtures at low 
temperatures. The theory is applied to the case of vibrating wires, immersed in the liquid. We 
present numerical results based on a direct solution of Landau-Boltzmann equation for the '^He 
quasiparticle distribution for the full scale of the quasiparticle mean-free-path I. The two-fluid 
nature of mixtures is taken into account in the theory, and the effect of Fermi liquid interactions 
and boundary conditions are studied in detail. The results are in fair quantitative agreement with 
experimental data. In particular, we can reproduce the anomalous decrease in inertia, observed 
in vibrating wire experiments reaching the ballistic limit. The essential effect of the experimental 
container and second-sound resonances is demonstrated. 

PACS numbers: 67.60.G-,67.10.Db 



^ I. INTRODUCTION 

o 

Vibrating wire resonators have extensively been used to study the properties of quantum fluids. Different regimes 
^ can be observed in Fermi systems at low temperatures. At relatively high temperatures the fluid obeys hydrodynamic 
description and the vibrating wire can be used to extract the density and the coefficient of viscosity of the liquid. As 
' ' the temperature is reduced below the Fermi temperature, the collisions between the quasiparticles, the elementary 
excitations of the liquid, get rare as the final states of a collision are limited by the Pauli principle. This leads to 
' increasing mean-free-path of the quasiparticles, and to increasing viscosity. With increasing viscosity the vibrating 
^ wire becomes more sensitive to its surroundings since more liquid is dragged to move with the wire, and this becomes 

'■^ limited by the experimental container. The increasing mean-free-path also leads to deviations from the hydrodynamic 
I behavior, which are first visible as a slip on the walls: the velocity of the fluid does not vanish at the wall but 
^ extrapolates to zero a slip-distance behind the wall (in the rest frame of the wall). Lowering the temperature further, 
^ the mean-free-path becomes comparable to the radius of the wire. In this regime the dynamics of the quasiparticle 
gas needs to be described by the Fermi liquid theory, which is similar as the theory of rarified gases but uses Fermi 

S distribution and has remaining interactions between quasiparticles. In the limit of lowest temperatures one enters the 
I ballistic limit where the mean-free-path exceeds the size of the experimental container. 

In this paper we provide theoretical calculations for a vibrating wire in the full range of mean-free-paths in a normal 
^ Fermi system. The theory is applicable to pure liquid ^He, but the range of applicability is limited by the superfluid 
O transition. Therefore the main application is in '^IIe-'*He mixtures, where the fermions remain in the normal state 
, ^ , in full range that has been measured. Several vibrating wire experiments have been made in ■^He-'*He mixtures. We 
compare our calculations to the measurements by Martikainen et a/.'-'-', which extend deepest into the ballistic regime. 
T-H In order to incorporate the '^He-component, we use Landau's Fermi- liquid theorjS generalized to the simultaneous 
I> presence of condensed bosons, as formulated by Khalatniko\3. A formulation of this theory adapted to the present 
OO work is given in Ref . A short account of the present calculations is given in Ref . [6] This letter highlights the Landau 
force, the macroscopic force caused by the Fermi liquid interactions. Here we calculate the Landau force among other 
effects and compare to experiments. We find strong effect of the experimental container on the resonance frequency 
of the vibrating wire. We also study sound resonances that affect the measurements at high frequencies or in large 
containers. 

•—^ Our calculation is the first two-dimensional solution of normal-state Fermi-liquid equations for general mean-free- 

I path. Even the cylindrically symmetric case we study, although effectively one-dimesional, is far more sophisticated 
• • than has been done before. Previous work has considered uniform bulk liquicP^^ or one-dimensi onal geometries 
• ^ such as an oscillating planar wall in semi-infinte liquicP^'^ or Poiseulle flow between parallel plane^^^*^. For more 
[v^ gen eral ge ometries calculations have been done only in limiting cases of small mean-free-patfP^El and in the ballistic 

^ limipim. 

We start in Sec.|ll]by stating more precisely the problem we study. We briefly present the equations we need, which 
are derived in detail in Ref. ^5, The symmetries of the problem are studied in Sec. |III[ In Sees. |IV| and |V] we study 
the special cases of vanishing mean- free-path in the hydrodynamic limit, and infinite mean- free-path in the ballistic 



limit. The numerical method is introduced in Sec. |VI| and the parameter values used are discussed in Sec. VII The 



results of the numerical calculations are presented and discussed in Sec. VIII 
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II. FORMULATION OF THE PROBLEM 



A. Main assumptions 

We study the force that a vibrating wire exerts on a Fermi hquid surrounding the wire. The wire is modeled as an 
infinite cyhnder. We assume that both the wire and the fiquid are in a stationary container that also has the shape of 
an infinite cylinder. These assumptions are not essential in principle, but they simplify the numerical calculations by 
allowing translation symmetry in the direction of the cylinder axis. The cross section of the wire is taken as a circle 
of radius a. The oscillation direction of the wire (chosen as the x axis) is perpendicular to wire axis (chosen as the z 
axis). 

For the container we consider two different shapes. By cylindrical container we mean a circular cylindrical container 
of radius b that is concentric with the wire. By slab container we mean that the wire is in the middle between two 
parallel plane walls spaced by 2h, and the oscillation direction of the wire is normal to the walls. Because of the latter 
case, circular symmetry is not valid in general. 

The force exerted on the liquid per unit length of the wire is denoted by F. We are interested in the linear response 
of the fluid. For that we express the force in terms of the mechanical impedance of the fluid, Z, as 

F = Zu, (1) 

where u is the velocity of the wire. In general Z is a tensor, but we concentrate on symmetric cases, where Z is 
a scalar. We assume harmonic time dependence cxp{—iujt) with angular frequency uj. Thus Z is complex valued, 
Z = Z' +iZ". The impedance is directly related to the resonance frequency /q and the line width A/ of the resonator 
by 

Z' = 27r\^p^Af, 

Z" = AiT^a^pUfo- fv.c). (2) 

Here pw is the density of the wire material and /vac the resonance frequency of the wire in vacuum. The restriction to 
linear response implies that the wire velocity is small compared to the Fermi velocity, u <^ vp, or any sound velocity 
in the fluid. Because of limitation to linear response, we systematically write all equations below in their linearized 
form. 

The plan is to calculate Z in the full range of the quasiparticle mean-free-path £, from the hydrodynamic regime 
£ ^ a to the ballistic limit £ ^ a. The liquid is assumed to be mixture of bosons and fermions, but the corresponding 
number densities and are arbitrary. Thus the theory applies to pure Fermi liquid as well (714 = 0). 

Besides linearity, another simplification is that we neglect the coupling between the normal and superfluid compo- 
nents. This can be justified in two different limits: at small frequencies and at small concentrations. In the former 
case we aim to calculate Z to accuracy that is correct to first order in auj /v-p . The coupling contributes to Z in second 
order in ui. Assuming that the relevant dimensionless parameter is auj/v-p, the coupling indeed is a small correction 
in the experimental case we study, where aw /v-p ^ 0.02 (Table |l]). In a finite container there appears a second length 
scale, h = 8a. Since hcu/vp is not very small, the small- frequency approximation is not necessarily accurate. For the 
superfluid component the low frequency approximation is well satisfled, huj/c <^ 1, due to the large sound velocity c 
of liquid ^He. We note that one more frequency dependent parameter, a'^ui/vp£, occurs in the theory, but it can have 



any magnitude (see Sec. IV I 



Alternatively, the normal-supcrfluid coupling can also be neglected in the limit of small "^Hc concentration. This 
is because the response of the superfluid to the motion of '^He is proportional to the ratio of their number densities, 
ns/n^ [S]. Thus, irrespective of frequency, the coupling can be neglected for small 113/ n^. 



B. Bulk Fermi-liquid theory 

The Fermi liquid theory was formulated by Landau to describe the low energy states of interacting Fermi liquidsP. 
It was generalized by Khalatnikov to simultaneous presence of condensed bosonS A new presentation of the Fermi- 
Bose-liquid theory is given in Rcf. 3 Here we briefly review the equations of Ref. [5] that are necessary for the 
numerical calculation. 

The central quantity in the theory is the quasiparticle distribution function ipp{r,t). It depends on the momentum 
p only through its direction p = p/p- In addition it depends on the location r and on time t. The distribution tpp is 
obtained from the more common quasiparticle distribution rip by integration over the magnitude of the momentum 
and making a transformation that partly decouples the normal and superfluid component&S-. The state of the bosons 
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is described by the deviation of the chemical potential S^4{r,t) = ^^{r,t) — ^ from its equilibrium value ^^^"^ and 
by the superfluid velocity Vs{r,t). 

One more important quantity is the quasiparticle energy shift on the Fermi surface, 5ep{r,t). It depends on 
and on ^/i4 as 



5/^4 + X! r 



{Flip ■ p'Hp')p' 



(3) 



21 + 1 



Here iiT is a parameter describing coupling to the superfluid component, Fi with ^ = 0, 1, 2, etc. are the Landau 
parameters describing interactions between quasiparticles. Pi are Legendre polynomials, and (. . .)p denotes average 
over the unit sphere of p. It follows from Eqs. ^ and (15) below that the coupling term contributes to Z proportional 
to oj^, and we neglect it. Since there are no experimental determinations of Fi for I > 2, we neglect the corresponding 
terms. Therefore Sep simplifies to 



Sep 



Fo 



Fi 



-p ■ {p''ipp')p'. 



l + Fo'^"'" ' l + Fi/3' 
The kinetic equation in the relaxation-time approximation takes the form 

:{ipp - Sep) + vpp ■ VV'p = --{ipp- ipp)- 



d_ 



(4) 



(5) 



Here v-p is the Fermi velocity that is related to the Fermi momentum p-p and to the effective mass m* by v-p = pp jm* 
In the collision term t ~ (./vp is the relaxation time and the local-equilibrium distribution 



V^J? = (^P'>P'+3p-(pVp')p' 



(6) 



Assuming time dependence exp(— zwt) and parameterizing the quasiparticle trajectories by r = ro + sp, the transport 
equation can be integrated in the form 



i^p{ro) ^ tppija + sop)e^''° + ds {ro + sp) ~ i — Sepir^ + sp) 

J so L« '"F 



(7) 



where k = 1/^ — iuj/vp. For convergence of the integral sq < 0, and the limit sq — > ~oo is approached. When a 
boundary is hit, the solution should be constructed piecewise, and the boundary conditions applied in between. 



C. Boundary Conditions 



The boundary condition appropriate for specular scattering from a surface moving with velocity u is 

^pp^2pp_2n{n-p)+2pp{n-p){n-u), (8) 

where n is the unit surface normal pointed to the liquid. In specular scattering, the incoming and reflected quasiparticle 
trajectories make the same angle to the surface normal on the wire. The boundary condition for diffusive scattering 
is that for outgoing quasiparticles (n • Pout > 0) 

2 

''Ppa^t = -2(r2. • Pi„-0p.„)p,„ + PF(Pout + • u, (9) 

where (. . .)p.^ is an average over half of the unit sphere (n • Pi^ < 0). In diffuse scattering the reflected (outgoing) 
quasiparticles on the surface of the wire are in equilibrium evaluated at the quasiparticle energy ep = vp{p — pp) + Sep. 
The distribution of reflected quasiparticles depends on the incoming quasiparticles only on the average. We can also 
consider mixed boundary conditions where fraction S of incoming quasiparticles is scattered specularly and fraction 
1 — S diffusely. We note that in using the linearized boundary conditions Q and Q, the displacement of the wire leads 
to a second order correction, which is neglected, and thus the boundary conditions can be applied at the equilibrium 
location of the wire surface. 

At the container wall, one possiblity is to assume diffuse scattering, i.e. to use Eq. ^ with u = 0. An alternative, 
simpler boundary condition is an absorbing wall. This reflects no quasiparticles and is described by the condition 
V'jJout = 0, where n • Pout > 0. Such a boundary condition could be a reasonable model for experimental cells that 
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have walls made of sintered silvei*-'-'. The absorbing boundary condition is also a theoretical tool to suppress sound 
resonances (Sec. VIII). 

We demonstrate the generalization of Eq. (j7| to include mixed boundary condition. By defining Xp{s) = ■(/'^(s)/^ — 
iujSep{s)/vp we write 

ijjpis = 0) = Sgce''''' + {2SpF{n^ ■ p)in^ ■ u) + (1 - S)[pf{p + ^n^) ■ u + g^]} e''"™ 

+ / dsXp{s)e^' + S ( " dsXp.{s)e'"' (10) 

for a trajectory that is reflected from the wire (surface normal n^u) at s = s^j, and hits the container (surface normal 
ric) at s = Sc- In the latter integral the prime in direction p' stands for the specularly reflected trajectory. We have 
defined two boundary condition terms 

gc{rc) = -2nc ■ (Pin'0pi„(r"c))pi„, (11) 

which depend on the locations and rc on the wire and container surfaces, respectively. For a trajectory with no 
collision with the wire, we get simply 

.0 

^p(O) = g^e'^"^ + / c^sXp(s)e'=^ (12) 



D. Force on the liquid 

The momentum flux tensor i^ 

n = p(o) Y + ^S^iJ + 3n3 {ppi>p)p, (13) 

7714 

where P^^^ is the equilibrium pressure. The superfluid density ps is defined as ps = min^ — Dm* / {1 + ^1/3), where 
1714 is the mass of a "^He atom, 77.4 the ^He number density, Z? = 1 — (1 + Fi/'ii)m^/m* and — pp/STr^h^ the number 

density of ^He. The force per area exerted by a surface element of the wire on the fluid is n • II evaluated at the wire 
surface — a{xcosd + ysm9) . Integrating this over the perimeter of the wire gives the force 



aps 
777.4 



/ Sp4,rid0 + 3an3 / {{h ■ p)p'ipp)p dO. (14) 
Jo Jq 

In linear theory both terms are proportional to it and thus this expression allows to determine Z ([T]). 

The two terms appearing in Eq. (14 1 can be interpreted as superfluid and normal fluid contributions, respectively, 
and correspondingly Z can be written as Z = Zg + Zn- We can calculate Zg as follows. Because the superfluid 
component is curl free, 'V x Vg — 0, the fluid flow can be described as potential flow, Vg — Vx. The ideal fluid 
equation of motioiP dvg/dt + VSp^/m^ = then gives 6114/7714 = —dx/dt. For small frequencies ui ^ c/a, v^/a we 
can neglect the compressibility of the superfluid and the coupling to the normal component, and assume V • t;s = 0, 
or V^x = 0. This has to be solved using the boundary conditions fi ■ Vg = fi ■ u on the wire surface and n -Vg = Q on 
the container walls. Thus the problem reduces to solving the Laplace equation. The solutions in several geometries 
have been found, see Ref. [H] The result is 

Zg ~ —ii^-Kc? pg G, (15) 

where the factor G depends on the geometry of the container. For infinite fiuid G = 1, for the cylindrical container 
G = (6^ + a2)/(fe2 ^ a^)^ and for the slab container G « 1 + -n'^a^ jYllr? 



In the limit £ ^ a the Fermi liquid theory reduces to the hydrodynamic theory (Sec. IV). In the extreme limit 
£/a ^ the viscosity is negligible, and the whole fiuid behaves like an ideal fluid except thin boundary layers on the 
walls. In this limit Z approaches 



^idcal 



= —iujna^pG, (16) 



where p = 771477,4 + 7^3773 is the total density of the fluid. It is convenient to present Z by its deviation from the ideal 
fluid behavior. Simultaneously, we define a dimcnsionless impedance Z by writing 

Z = 2'idoal + a773PF^- (17) 
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Using the symmetry assumption || m on equation (14) and comparing to (17) gives 

^^T^W7^— + / {in-p)iP-u)i:j,)j,d0, 18 

1 + Fi/3 VF pfu Jo 

where the second term remains to be calculated numerically. It follows from the equations above that Z depends on 
the dimensionless parameters auj/vp^l + Fi/3), £/a, S, Fq, Fi, and b/a in the cylindrical container or h/a in the slab 
container. Note that choosing + Fi/3) as an independent variable instead of 17 = auj/v-p is convenient because 
this combination frequently occurs in the hydrodynamic region, an example being the first term in Eq. ( |18[ ). 

III. SYMMETRY 

The problem stated in Sec. |II A| has symmetries that can be used to simplify the numerical calculation. We need to 
consider scalar and vector functions of location r. Because of the translational symmetry, there is no dependence on 
z, and therefore effectively r = xx + yy, or using polar coordinates r = r{x cos + y sin 6). Because we are studying 
linear response, these functions depend linearly on the wire velocity u. In the slab container we have reflection 
symmetries in the planes perpendicular to both symmetry directions x and y. Choosing u ~ x we have for scalar 
functions / and vector functions / — frr + fgO the symmetries 

f{r,e) = /(r,-^^)--/(r,^-0)--/(r,0 + 7r), 
U{r,e) = U{r,-e) = -}r{r,7:-6)^-U{r,e + T:), 

fe{r, 9) = ~fg{r, -9) = fg{r, n - 9) ^ -fg{r, 9 + n). (19) 

For the cylindrical container we argue as follows. Because of rotational symmetry in simultaneous rotations of 
r and u, a scalar function has the form f{r,u) = f{r,u,9), where 9 denotes the relative plane angle and r = \r\ 
and u ~ \u\. Assuming reflection symmetry with respect to the a;-axis f(r,u,9) — f{r,u,—9). This allows to write 
/(r, f ■ u,u) or f{r, r ■ u,u). Requiring / to be linear in u limits this to the form 

f{r,u) = f{r)r-u (20) 



with some function /(r). Note that the symmetry (20 1 often appears in solutions of the Laplace equation. Our 
derivation shows that this symmetry remains valid in the case of an arbitrary linear and isotropic equation. 

Next, we consider a vector field f{r,u) — ffr{r,u) + 9fg{r,u). Rotational symmetry limits the two functions to 
the forms fr(r,u,9) and fg(r,u,9). Assuming reflection symmetry fr{r,u,9) = fr(r,Uj —9) allows to write this into 
the form fr{r, f ■ u, u). Assuming reflection symmetry /e(r, u, —9) — fg{r, u,tt + 9) allows to write this into the form 
fg(r,6 ■ u,u). Finally, assuming linearity with respect to u implies 

fir,u)=fr-ufrir) + ee-ufe{r). (21) 
The results above are applied to the bulk angular averages 

c(r) - (^p(r))p, b{r) = 3{p^Pp{r))p (22) 



and to the boundary-condition averages gw{fw) and gd'i'c) (Hi- The functions c and b are needed to calculate Scp 
lil and 4,^^ ([6 



IV. HYDRODYNAMIC LIMIT 

The Fermi-Bose liquid theory reduces to two-fluid hydrodynamic theorjEH in the limit of short mean- free-patlPSl. 
The normal fluid component is described by a Navier-Stokes equation where the density is the normal fluid density 
Pn = m*n3/(l -|- i^i/3) = p — ps and the coefhcient of viscosity r/ = ^n^p-pL The diffusive boundary condition leads to 
no-slip boundary condition, which means that the fluid velocity at a wall equals the velocity of the wall. The specular 
boundary condition leads to perfect slip, where the transfer of transverse momentum between the liquid and the wall 
vanishes. Applied to the wire surface this means Tire = 0. The hydrodynamic equations have analytic solution in 
some cases, which we discuss below. 

In the low-frequency limit the normal and superfluid equations decouple and both components can be considered 



as incompressible. In this limit, the superfluid component Zg (15) was found above. Analytic solutions for the normal 
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FIG. 1: The impedance Z = Z' -\-iZ" in the region of short mean-free-paths, £ < a. The curves are plotted using £ as parameter. 
The hydrodynamic results in unlimited fluid for diffusive wire ( 23 1 and for specular wire (|24[) are shown by thin dashed and 
dotted lines, respectively. The slip correction to the former is shown by dash-dotted lin#^. The hydrodynamic result for a 
diffuse wire in diffuse container of radius 6 = 8a is shown by thick dashed line. Also shown are our numerical results: a large 
cylindrical container with a diffusive wire (green) and specular wire (pink) with absorbing walls, a small container, b — 8a, 
with a diffusive wire (light blue) and specular wire (light green) in diffusive container, and a slab of thickness 16a (blue) with 
diffusive walls. The parameters are appropriate for a 3.6% mixture, aoj/vp — 0.0180, Fq — —0.25, and Fi = 0.22. Experimental 
data from Ref. I is shown for reference (yellow points). 



component are known in the following cases. The no-slip case in unlimited liquid was solved by Stoked^. The result 



a? H^^^ (qa) 



(23) 



Here, H^''\x) are the Hankel functions with complex argument, prime indicates derivative, q = {I + i)/S, and 
S = y/2ri/pnUj is the viscous penetration depth. The no-slip case in a cylindrical container was calculated by Carless, 
Hall, and HoolPl. The case of perfect slip in unlimited liquid was calculated by Bowley and Owers-BradlejP^, and 
the result is 



-iujTTa^ Pn[l + 



2aqH^o^\qa) - a^q^H'o'' (qa) 



r(l) 



(24) 



The hydrodyn ami c resu lts a re illustrated in Fig. [T] 

The results (23) and (24) illustrate the argument made in Sec. 
,2 /r2 _ 572, 



HA 



that Z (171 is linear in auj/v-p but one should 



allow /5'^ ~ ^a'^uj/vF£{l + |i^i) to have arbitrary magnitude. In the parameter space formed by auj/vp and i/a, Z 
is singular at the point (0,0). 

It is possible to go beyond the low-frequency approximation at the expense of neglecting viscosity. We study the 
case of the cylindrical container. At frequency w the velocity potential x of ideal fluid should satisfy the wave equation, 



where c is the sound velocity. We make the ansatz 

X{r,0) - [AMkr) 



BYiikr)] cos 61, 



(25) 



(26) 



where Ji and Yi are Bessel functions and here k = lo/c. The radial velocity has to equal u cos 9 at r — a and vanish at 
r — b. This gives conditions from which A and B can be solved. Calculating the force similarly as for the superfluid 
component in Eq. (14), we get 



deal 



iirapu Ji{ka)Y{{kb) - Yi{ka)J[{kb) 
k J[{ka)Y{{kb) ~ Y[(ka)J[{kb) ' 



(27) 



where p is the density of the fluid. To first order in lo this reduces to the low- frequency result Zjdoai (16 1. Sound 
resonances are found at frequencies where the denominator of ^idoai (27) vanishes. For large bja this condition reduces 
to zeros of J[{kb). The lowest zero appears at kb = 1.84. 
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The two-fluid hydrodynamics allows two sound modes. The first sound is in-phase motion of the normal and 
superfluid components. The corresponding sound velocity is close to the sound velocity of pure liquid ^He, which 
is relatively high, and is not of interest here. Second sound is the counter motion of the normal and superfluid 
components. In "^He-^He mixture it is essentially a compressional mode of the normal fluid as the superfluid motion 
cause s onl y a small correction to the sound velocity because of smallness of 713/714 In the present approximation 



(Sec. II A) the sound velocity is vpy^{l + Fi/3)(1 + Fq)/3 w OAvp in the absence of dissipation. This result will be 
applied below to identify resonances of second sound. 

An extension of the hydrodynamic theory to slightly longer mean-free-path is known as slip theorj033. The lead- 
ing corrections appear in the boundary conditions. Considering the diffusive boundary condition (|9]), the outgoing 
(scattered from the wall) quasiparticles are in equilibrium with the wall, but the incoming quasiparticles generally 
are not. In the case of short £ these populations get rapidly mixed and therefore also the incoming quasiparticles are 
in equilibrium with the wall, which leads to the no-slip boundary condition. With increasing i this is no more the 
case. In slip theory one uses bulk hydrodynamic theory together with boundary conditions that assume the velocity 
of the fluid to extrapolate to the wall velocity the slip length C behind the wall. Microscopic calculation of ^ for 
degenerate fermions on a planar wall was made by Jensen et al!^. The generalization of the slip boundary condition 
to curved surfaces and partial specularity is discussed in Refs. I16|17|24l The slip increases the validity range of the 
hydrodynamic theory to slightly larger £, see Fig. [T] Several extrapolations of the slip theory for vibrating wires have 
been suggested, see Carless, Hall and HoolP^, Guenault et al^, and more recently by Bowley and Owers-Bradlejff^ 
together with Perisanu and Vermeulerfi^. None of these have attempted to include the effect of the container or the 
Fermi-liquid interactions. 

V. BALLISTIC LIMIT 

At the lowest temperatures one enters the ballistic regime, where the mean-free-path of quasiparticles becomes 
large compared to the experimental dimensions. In the extreme case one can neglect the collision term in the kinetic 
equation ([s]). In this ballistic limit the quasiparticles still interact thro ugh the Fermi-liquid interactions, which appear 
through Sep. In spite of this, some results can be obtained analyticalljffSMEol^ 

If, in addition to the collision term, one also neglects the time dependence, the kinetic equation ([5| reduces to 
p ■ V'(/'p(t) = 0. This implies that ipp is constant along trajectories, and changes only when the trajectory hits a 
wall. If the fluid far from the wire is in equilibrium, one finds V'p = on incoming trajectories. The distribution 
for outgoing trajectories is then obtained from the boundary conditions ([s]) or ([9]). The impedance is obtained by 
integrating over the wire surface, and it giveJ^^El 

^diff = -^anspF, 

^spcc = -^arispF (28) 

corresponding to the diffusive and specular boundary conditions. These are purely dissipative since our assumption 
of w — > 0. These results have been generalized to the cylindrical container by taking into account the reflection of 
quasiparticles from the container wall back to the wire in Ref. [20!, 

At finite frequency the main complication arises from the interaction term Sep in the kinetic equation ([s]) . Neglecting 
this term it is possible to solve the ballistic limit at arbitrary frequency in the cylindrical container. We assume diffusive 



boundary condition on both the wire and the container. Using definitions (11) and the symmetry (20), the boundary 
conditions ^ can be written as 

2 

V'jjw.out (."^-w) = (Jw-n ■ u + |5F(Pw,out + 7^n) ■ u, (29) 
V'pc.out('^c) = gcfic ■ u, (30) 
where and gc are constants. Tracing the different trajectorices one can now calculate the averages appearing in 



definitions (11). We get g^, = Dijc and 



2 2 

gc = Age + B{g^ + -pp) + -Cpp- (31) 
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The coefficients A, B, C, and D are expressed as integrals 



A = 



B = 



dCsin^C 







arcsin(a/6) 



(^7 cos 7(1 — 2cos^ 7) exp 



/ 2iujb cos 7 
\ vp sin ( 



a 2 
b 



exp 



iuj{\/b'^ — a^x^ — a\J\ — x^) 



sin 



C 



D 



y— [ dCsin'^^ / dxW 1 — — a;2 exp 
b-n- Jo Jo V 



j{^/b'^ — d?-x'^ — a\J\ — x^) 



tip sin C 



d^sin C / c?x -a; 
Jo \b 



62 



a;^ exp 



iuj{^/b'^ — a^x"^ — aVl — x^) 
sin ( 



where C is the angle between the trajectory and the cylinder axis. We can now solve (31 1 and get 



|(g + C)pp 
1-A-BD' 



Continuation to calculate the force gives 



7ran3 



43pF 
48 



1 



2a 




36 1 - A 



(32) 



(33) 



(34) 



In the special case w 



the integrals ( 32 1 reduce to 
1 




(35) 



where D still needs to be calculated numerically These reproduce the result given in Ref. 

The result ( 34 ) is illustrated in Fig. [2j We see that Z„ oscillates around its static value approximately like the 
exponential oc exp(2ia;6/wp). The reason is that the quasiparticles excited by the wire are reflected from the container 
wall back to the wire, but delayed by time ~ 2b/v-p. The first constructive interference (after the zeroth one at = 0) 
corresponds to w « Trvp/b. This is approximately a factor 4 larger than the first second-sound resonance frequency 
in the hydrodynamic limit (Sec. IV I. The difference arises because the second sound velocity is smaller than the 
quasiparticle velocity vp, and because the first zero of J[{kb) is at kb = 1.84, which is less than tt. Second sound is 
strongly damped in the ballistic regime because of excitation of quasiparticles, known as Landau damping. 

The results above could be generalized to include the Fermi-liquid interactions. However, the integrals get very 
cumbersome even in the case of unlimited liquid. Instead of partial analytic solution, we therefore prefer the full 
numerical solution in the following. 



VI. NUMERICAL METHOD 



We solve the quasiparticle distribution i}}p{r) in a discrete grid around the wire and for discrete momentum directions 
p. The location r„ = xrcoad + yrsin^ can be parameterized with cylindrical coordinates r and 0, where we have 
fixed X — ii. The momentum direction pji = sin (^j(cos/3/£ -I- sin/3;y) -I- coaQz is parameterized by angles C,j and Pi. 
For the cylindrical container we use a grid with exponentially increasing spacing Ar in the radial direction, and a 
fixed spacing A0. In the slab container, we use a cylindrical lattice near the wire, and rectangular lattice near the 
container walls, with partial overlap between the lattices. In y-direction the length of the slab is increased until the 
results converge (absorbing walls are used). 

The quasiparticle distribution il^nji = '^Ppjii'^'n) for each lattice point n and for each discrete trajectory direction pji 



can be written, using Eq. (10 1, as the sum 



l-njl 



Bnjl + S Hnjl + (1 — S)Qnjl, 



(36) 
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FIG. 2: The normal-fluid impedance Zn = Z' -\- iZ'^ in the baUistic hmit in the cyhndrical container of radius b = 10a. The 
impedance Zn is plotted relative to its static value Zq = (437r/48)an3PF (28 1 in unlimited fluid, and the total impedance Z is 
obtained by adding Zs ( 15 1. The solid lines are full numerical result, the dashed lines the approximation ( 34 1, and the difference 
between them is the Landau forc^. (a) The real (red) and imaginary (blue) parts of impedance as a function oi Q, = aLo / v-g . 
(b) Parametric plot of Zn in complex plane with f2 as a parameter. In both panels diffusive boundary conditions are used and 
the fuU calculation uses Fo = —0.44 and Fi = 0.449. 



where /„ji = /° ds Xp^., (s)e'''* + 5* /j*™ ds Xp^.^ (s)e'^'* is the mtegral part, i?„ji = [2,SpF{n-p){h-u)^-{\ — S)p-p{p+'^n) 

Injl 



u\e^^™ is the i/j-independcnt part of the boundary condition, Hnu describes the ^-dependent part of the boundary 



condition on the container, and Qnji describes the boundary condition on the wire. The integral is calculated from 
r„ along the trajectory pji, \sw\ is the distance to the wire (if there is a collision with the wire at point r^), and \sc\ 
is the distance at which the container wall is hit, at point re- If there is no collision with the wire, we get from Eq. 
( [T2| ) i}}n]i = Inji + Hnji, whcrc Inji = ds Xp.^{s)e'"' . We write gc.t = gc{rt) and g^^t 9w{rt) at each lattice point 



rt on the surfaces, and define 

t 

(37) 

t 

where the sums are over all lattice points on the corresponding surface, H^tji — Wtc'"'^ and Qntji = Wte^^^ . The 
weight factors Wt are obtained by interpolating from the lattice points closest to or on the surfaces. In the 



following the averages (22) are denoted by '^n^ — Cn, ^'1^'' — 6r, n, and 'i'n^ = hg^n- Now, the integral part Inji can 



be conveniently presented as the sum 

^n.^ = EE«itU*L^^ (38) 

h=l m 



where the coefficients a^^mji weights with which the equilibrium distribution at other lattice points r„i affects 

ipnji- The coefficients are calculated by numerical integration, using a modified Simpson's rule, which takes the 
exponential term exp(fcs) into account exactly. It is necessary to use interpolation between lattice points. Taking the 
numerical average over the trajectories, we can write the averages ( |22[ ) as a sum 

= E + + E ^n* ff^.* + E (39) 

m,/i t t 

where 

^ih _ I Ah) , 

df^ = Se^-{p,ialii^i)p,,, (40) 
and similar averaging gives Bn^ , h!^^ , and qI^I . The diffusive boundary condition terms can be written as sums over 
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1.8% 


3.6% 


5.6% 


6.6% 


7.0% 


9.5% 




1.8% 


3.6% 


5.6% 


6.6% 


7.0% 


9.5% 


m* jm-i 


2.26 


2.31 


2.34 


2.35 


2.60 


2.63 




2.45 


2.50 


2.53 


2.54 


2.85 


2.88 


iiF(m/s) 


21.10 


25.97 


29.60 


31.07 


29.56 


32.36 




19.46 


24.00 


27.38 


28.75 


26.97 


29.55 


anapp (kg/ms) 


0.0058 


0.0145 


0.0259 


0.0321 


0.0394 


0.0589 




0.0058 


0.0145 


0.0259 


0.0321 


0.0394 


0.0589 




0.151 


0.219 


0.266 


0.284 


0.269 


0.301 




0.138 


0.201 


0.245 


0.261 


0.243 


0.272 


auj/vY, wire 1 


0.0222 


0.0180 


0.0158 


0.0151 


0.0159 


0.0145 




0.0241 


0.0195 


0.0171 


0.0163 


0.0174 


0.0159 


aLj/vF, wire 2 








0.0233 


0.0245 


0.0224 










0.0252 


0.0268 


0.0245 



TABLE I: Parameter values used in calculations at different concentrations X3. The first four concentrations are at saturated 
vapor pressure, the last two at P = 10 atm. We use a = 0.284 at saturated vapor pressure and a = 0.207 at P = 10 atrrP^ 
The frequency parameter fl = au/v-p is calculated for vacuum frequency / = 1202.85 Hz (wire 1) and / = 1857.7 Hz (wire 2). 
On the left we have used the hydrodynamic mass vriH/'m^ — 2.15 as a starting point, and on the right mH/ms = 2.34 (see 
text). 



lattice points m as well, 



9c,t = J2 fSm*™^ + Bc,t + H,^ts9c,s + J2 Qc,ts9^.s 
m,h s s 

9w,t = gwlrn^rn + H^^tsQc^s, 



(41) 



where the new coefficients g^'trn' ^cti H^^tsi Qc.ts, 9wirm ^^'^ H^ ts are similar to the terms above, but the averages 
are taken over the incoming trajectories only. Once the coefficients have been calculated, we can solve the averages 
( 22 1 from the matrix equation 



m.h 



* = i:>* + B^* = (/ 



D)-^B. 



(42) 



If the total number of lattice points is iV, the number of points on container surface is V , and the number of points on 



wire surface is T, then \1/ is a 3iV - 



- y + T component vector 4- = (^'i'^^'^^^ . . . • • ■ , 9c,v ■, gw,i, ■ ■ ■,gw,TY\ 

, Bc,v,0, 0)^, and the (37V + V + T) x (37V + V + T) matrix D contains the 



B = {B[ 
terms dl'i 



nm^ nt ' 



r(i) 



M Ah) 
' nt ' Hc^tm 



H, 



(h) 



c,tsi '^c,tsj 9w,tm' 



and Hyj fg in appropriate order. 



When the averages are known, we can calculate the full distribution ijjnii on the wire surface using Eq. ( 36 ) . The 



force exerted by the wire on the fluid is then readily calculated from Eq. (18 1. If cylindrical symmetry is assumed, the 
treatment is similar, but since the dependence of the averages on is of a simple form [(20 1 and (21 1], the calculation 
becomes effectively one-dimensional. 



VII. PARAMETER VALUES 



In order to relate the numerical calculations to experiments we need to consider the parameters that characterize 
•^He-^He mixtures and the experimental setup. 

The molar volume of a mixture of molar concentration 2:3 — N^/ {N^ + N4) at temperature T and pressure P i^^H 

V^{x3,P,T)^V4iP,T)[l + a{x3,P,T)x3], (43) 

where V4(P, T) is the molar volume of pure '*IIe and a is the so-called BBP parameteJ^^. For the molar volume of 
pure ■*He, we use the data of Tanaka et al^, and for a{P) at zero temperature we use the results of Watson et alW^. 
The number density of the ^He component is obtained from 

"^=Vi(l-Ka.3)' ^''^ 
where a is the Avogadro number. From ^3 we can calculate the Fermi momentum 

PF = ^3^2^, (45) 

If we know the ^He quasiparticle effective mass m*, we can also calculate the Fermi velocity v-p = p-p/m* , which 
is needed in calculating the dimensionless frequency parameter = auj/v-p. The effective mass m* is calculated 
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1.8% 


3.6% 


5.6% 


6.6% 


7.0% 


9.5% 


Fo, Ref. |36|37| 


-0.17 


-0.25 


-0.33 


-0.36 


-0.38 


-0.47 


Fo, Ref. EH] 


-0.12 


-0.21 


-0.28 


-0.32 


-0.24 


-0.29 



TABLE II: Fo at different concentrat ions b ased on data presented in Ref. [35]. The first fine gives results from a fit to the 
osmotic pressure data of Landau et alW^^, and the second line gives a fit to the second sound velocity data by Murdock and 
CorrucinP^. The first four concentrations are at saturated vapor pressure, the last two at P = 10 atm. 



using the interpolation formula by Krotscheck et al P^^^^ [ As a starting point, the zero concentration limit effective 
mass (so called hydrodynamic mass) mn is needed. We have used two alternative values: mn /m^ — 2.15 based on 
extrapolation of the data by Simons et al.^^ following Krotscheck et alfi^, and rriH/m^ — 2.34221124] ^ero pressure. 
Here 7713 is the mass of a '^He atom. At higher pressure P = 10 atm we have used mn /m^ = 2.39 and mH/m^ = 2.64, 
correspondingly. Numerical values are given in Table |l) For Fi we use the approximative formula Fi = 3{m* /niH — ^)- 
For Fq we use a fit of the form Fq = Axl^^ + BxV'^ to two sets of data presented in Ref. |35l see Table II 



3 r ^-^^ 

4/3 .x,riT 



We present Z (11] using the dimensionless Z in order to remove the strong dependence oc rig on the He density. 
The experimental result J^l^ are given as (/o. A/) pairs. In order to transform them to Z we use Eqs. ([2]) and ([17]), or 



Z' = ?![!f^A/, 
ansPF 

Z" = ^^^^[(l+jG^)/o-/vac]. (46) 

Several vibrating wire measurements have been done in "^He-^He mixture j^^ * ^'^ * ^^ * '^^ . Here we study the data of 
Martikainen et al\^, which extends deepest into the ballistic regime. The parameters needed for the conversion 



(46) are the wire radius a = 62 /im, the distance between the container walls 2h = 16a, the vacuum frequencies 
/vac = 1202.85 Hz for wire 1 and /vac — 1857.7 Hz for wire 2, the density of the wire pw — 16700 kg/m'^ and the 
concentration-dependent densities as explained above. The vacuum frequencies /vac were measured before filling the 
cell. The measurements with filled cell show small deviation from the predicted high temperature behavior as if the 
vacuum frequency had some variation during the measurements. We have compensated this by adjusting /vac for each 
concentration separately in the range ±0.2 Hz so that at high temperatures Z extrapolates to zero. 



VIII. RESULTS AND DISCUSSION 



In this section we present the results of numerical calculations for the full scale of the mean-free-path In order 
to show the dependence of Z on different parameters, we select a "basic set" of parameters: ^2/(1 -|- fi/3) = 0.0145 
(where 51 = auj/vp), the cylindrical container with b = 10a, S — 0, Fq — —0.33, Fi = 0.266, and diffusely scattering 
container walls. Then we vary each of the parameters + Fi/3), b/a, S, Fq, and Fi separately while keeping the 
others fixed at the basic set values. The real and imaginary parts of the impedance Z are interpreted as dissipation 
and shift of resonance frequency, see Eq. ([2]). 

The dependence of Z on 51 = auj/vp is shown in Fig. |3] In the chosen dimensionless variables, the dependence on 
frequency uj appears through this parameter only. The reactive part Z" vanishes for vanishing frequency. For the 
smalle st va lues of Jl the dependence seems nearly linear, but the hydrodynamic limit is more complicated as discussed 
in Sec. IV The dissipative part Z' has weaker dependence on frequency and reduces to a finite value at Jl — >■ 0, except 
in the limit £ — ^ 0. At larger ft nonlinearity appears in Z" at all £/a. The last curve at Jl = 0.08 differs essentially from 
the others. This behavior is caused by buj/vF approaching unity, where the lowest second-sound resonance appears in 
the container. More precisely, the lowest pole of Eq. (27 1 gives the resonance at flc = 0.0890. The impedance in the 
neighborhood of the resonance is shown in Fig. |4] Resonance features that resemble our results have been observed 
in experiments using quartz tuning forks, that work at considerably higher frequencies (33 kHz) [5]. 

The resonances in Fig. [l] appear strongest at small £ while they are damped at larger £. The behavior of the 
ballistic limit point is shown in Fig. [2j Instead of sharp resonances one sees oscillatory behavior. Thus we see 
complete change-over from the second sound resonances in the hydrodynamic regime to quasiparticle interferences in 
the ballistic regime. 

The effect of confinement is studied in Fig. [5] in terms of the radius b of the cylindrical container. For diffusive 
container walls b has a strong effect on the results: the dissipation and the minimum frequency are larger for smaller 
containers. If we use absorbing walls for the container, the container size does not have such a drastic effect. It is 
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FIG. 3: The effect of the frequency parameter Q — olo/vf. In this and the following figures 4|8 panel (a) shows Z' vs. Z" 
with i/a as a parameter in the curves. Panels (b) and (c) show Z' and Z" as a function of £/a, the mean-free-path scaled by 
the radius of the wire. The variation of fl can be understood as change of concentration, or change of the vacuum frequency 
of the wire. The dissipative part Z' has a constant part independent of Q, but the reactive part Z" vanishes for — >■ 0. For 
other parameters, the basic set is used; b — 10a, S = 0, Fq = —0.33, Fi = 0.266, and diffusive container walls. The curve for 
the largest value fl = 0.08 differs essentially from the others because it is close to the lowest second-sound resonance, see Fig. [4] 
for more details. Here and in Figs, 
values of Q, and Fi. 



4)8 the dashed black curve gives the hydrodynamic (HD) Stokes solution ( 23 1 for basic-set 
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FIG. 4: Near the expected second sound resonance fic = 0.0890 we observe peculiar behavior in both dissipative and reactive 
parts of the impedance. The three short-dashed lines with Q — 0.06, 0.07, and 0.08 show the usual behavior below the resonance. 
The three solid lines, Q = 0.085, 0.088, and 0.089 are close to the resonance and make loops in the Z' — Z" plane as seen 
in panel (a). The long-dashed curves, — 0.090, 0.091, and 0.095 show the behavior above the resonance; it appears that 
these curves would form loops as well, if we could extend the calculations to smaller I. For the resonance curve, Qc ~ 0.089, 
dissipative part has a maximum at around £ = 0.007a and a minimum at ^ ~ 0.08a, while the reactive part has a maximum at 
£ ~ 0.015a and a minimum at £ ~ 1.2a. These extremal points move towards smaller £ with increasing Q. 



understood that the reflecting waUs aUow for second-sound resonances, which occur at bfl/vp ^ 1, and have a large 
efi^ect at higher frequencies or at larger container sizes. 

The dependence on the specularity parameter S of the wire surface is shown in Fig.|6] Increasing S means increasing 
slippage on the wire surface, and leads to less fluid moving with the wire. Then the dissipation and the change of 
resonant frequency are smaller than for fully diffuse wire. We emphasize that the parameter S is assumed to be 
independent of £. We see that the dependence on S is nearly linear except in the hydrodynamic region, where the 
fully specular case 5 = 1 stands out from all other values S < 1. For absorbing container walls the effect on the 
resonant frequency in the ballistic limit is small. 

The contribution of the quasiparticle interactions to Z is called the Landau forc^. The dependence on the inter- 
action parameter Fq is studied in Fig. [7j The effect of Fq appears via the change of density. In the hydrodynamic 
regime the density stays almost constant and therefore the effect of Fq is small. In the ballistic regime density changes 
are essential and Fq has a large effect on Z" . For stability Fq has to be larger than —1. The effect of Fq on Z" seems 
to be amplified for Fq approaching the instability. For a low frequency such as in Fig. [7| the effect of Fq on the real 
part Z' is small. With increasing frequency the main effect turns from Z" to Z' , as can be seen in the ballitic- limit 
Fig. [2] In the large frequency limit the main effect of the Landau force is the reduction of Z' (for Fq < 0). 

The effect of the second interactions parameter Fi is shown in Fig. [sj For constant fl/ (1 + Fi/3) its effect is very 
similar to that of Fq, the curves almost coincide. The difference is that increasing Fi is compensated by a smaller 
decrease of Fq. The opposite tendency can be understood based on the qualitative explanation of the Landau force 
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FIG. 5: The effect of the container radius h. The solid lines correspond to diffuse container walls, while the dashed lines 
correspond to absorbing walls. For small I the container radius has only small effect, and the curves merge. At larger the 
maximum shift of frequency is larger for larger container, while in the ballistic limit the maximum dissipation is larger for 
smaller container. The largest radius corresponds to effectively unlimited fluid. We see that although Z" > for h = 100a in 
the ballistic limit, it is quite small in comparison to the smaller containers. In the case of absorbing walls, the curves end to 
the same point in the ballistic limit, close to the end point of the large container curve with diffuse walls. 
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FIG. 6: The effect of the specularity S of the boundary condition at the wire surface. For small i the effect of S is hardly visible, 
except for the purely specular curve 5* = 1 which stands out in the hydrodynamic limit for Z" . At larger t the dependence on 
S is uniform in both Z' and Z" . Both the maximum frequency shift and the maximum dissipation are larger for a diffusive 
wire. 
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FIG. 7: The effect of Fq. In the Z' vs. Z" curve, -Fo controls the end point and the slope of the ballistic branch of the curve, 
while its effect is small in the hydrodynamic region. From panel (b) we see that the effect on the dissipation is small, apart 
from the smallest value shown here, = —0.8 in the ballistic regime. The effect on the frequency is larger, as discussed in 
Ref. [6l The vertical (a) and horizontal (c) dashed lines mark the ballistic limit value of Z" in the case of Fo = 0. 



given in Ref. [S) Namely, the oscillating wire creates a beam of quasiparticles that affects the quasiparticles that are 
incident on the wire. The relevant interaction thus has p ■ {p' < in Eq. Q and thus Fq and Fi appear with weights 
of opposite signs. 

We notice that the dissipative part Z' depends weakly on Fq, Fi, and fi/(l + -^1/3) (except close to resonance), 
while the reactive part Z" has a stronger dependence on the parameters. For S the both parts are affected in a similar 
fashion, and for b the dependence is rather complicated. 

Our calculations are compared to experiments of Martikainen et al\^ in Fig. [o] We model the experimental 
chamber as the slab container with diffusive boundary condition on both the wire and the container. Since there is 
uncertainty in m* and Fq, we show four combinations of parameters for each concentration. We also show results for 
the absorbing boundary condition at the container walls and for the cylindrical container of radius b — 8a with the 
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diffusive boundary condition. Although the agreement between the experiments and the calculations is not perfect, 
it is for most cases satisfactory, noting that no fitting parameters have been used. 

In the ballistic limit the resonance frequency increases to value higher than in the high temperature limit. This 
overshoot was analyzed in detail in Ref. HI There the overshoot was divided into four different contributions: 1) The 
■^He part of the fluid decouples from the ideal fluid flow around the wire. 2) Part of "'He moves with the quasiparticles 
and therefore is also decoupled. [These two contributions correspond to the normal density p„ being subtracted from 
the total density p appearing in (jlS])]. 3) The Landau force due to quasiparticles interactions, caused mainly by Fq^ 
adds elasticity to the Fermi liquid and thus increases the resonant frequency. 4) There are corrections caused by the 
finite size of the container, in particular the effect of quasiparticles reflected from the container back to the wire ( 34 1 . 
The four contributions are shown by vertical lines in Fig. [9j 

We see that for diffuse chamber walls the calculated overshoot is larger than in the experiments. It is possible that 
the ballistic limit was not quite reached in the experiments. Alternatively, the calculations with absorbing container 
walls give overshoot of the same magnitude as the experiments. The experimental cell was surrounded by porous 
sintered silver, which is likely to absorb some of the quasiparticles rather than reflecting them. We note that the slope 
of the ballistic branches of the calculated curves would better fit to their experimental counterparts if larger values of 
—Fq were used. Also, allowing non-zero specular-scattering fraction S would scale down the calculated curves, giving 
better fit for some of the concentrations, but a systematic improvement is hard to obtain. We point out that for the 
calculations in cylindrical geometry we have fixed the radius b = 8a rather arbitrarily: using 6 as a fitting parameter 
would certainly lead to better agreement with experiments, but again it is difficult to find a single value that would 
fit all the concentrations. 

We see from table [l] that for concentrations 5.6% and 7.0% for wire 1, the frequency parameters are nearly the 
same, 57 w 0.016 (or ft w 0.017 depending on the choice of mn)- However, the experimental results, converted to 
the {Z" , Z')-curve, differ considerably from each other, as seen in Fig. 10 1. The corresponding numerically calculated 
curves differ only slightly, due to small difference in Fq. There seems to be some problem involved with the two curves 
at higher pressure for wire 1, and for these the fit to numerical calculations is by far poorest of the nine cases, as seen 
in Fig. and f. 

Another interesting coincidence is that for the 1.8% curve of wire 1 and for the 9.5% curve for wire 2, the frequency 



parameters are similar 17 « 0.022 (or 17 0.024). The experimental curves shown in Fig. 10 3 are nearly identical 
in the hydrodynamic region, and only start to differ at £ > a. The slopes of the ballistic branches of the curves are 
different, which is conveniently explained by the different values of Fq for the two concentrations. 

There is a systematic difference between the results of wires 1 and 2. For wire 1 the calculated minimum frequencies 
are smaller than experimentally observed, while for wire 2 the opposite is true. As already discussed, the high pressure 
curves (7% and 9.5%) for wire 1 seem anomalous. There also seems to be fine structure in the measured curves that 
is not present in the calculations. In particular, the nearly straight part from the bend of the curves to the ballistic 
limit seems to have small positive curvature in the measured curves. As shown above, the theory gives second-sound 
resonances, but at higher frequencies or at larger wall distances than used in the measurements. In other directions 
the experimental container is ten times larger and therefore second-sound resonances in these directions should take 
place at the experimental frequencies. Ideally, such modes are not coupled to the oscillation of the wire, but there 
might be coupling if the wire is not perfectly aligned with the slab. This may be the origin of the observed structures, 
which are not reproduced by our ideal infinite-cylinder model. 

Our calculation can be generalized in a couple of relatively simple ways. We can tilt the oscillation direction in the 
slab, we can allow cylindrical containers of more complicated cross section, and we can allow two relaxation times 
(one for / = 2 spherical harmonics of ipp and another for higher harmonics). We have tested all these, but they do 
not seem to give any obvious improvement in the comparison above. Further possible generalizations could be more 
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FIG. 9: Comparison of theory to measurement^^. Wire 1 has been measured at six concentrations (1.7%, 3.6%, 5.6%, 6.6%, 
7.0% and 9.5%) and wire 2 at the three highest of these. At concentrations 7.0% and 9.5% the pressure is 10 atm, while at the 
others it is saturated vapor pressure. The experimental results are shown by blue data points. The principal theoretical result 
is the red solid line in each panel. It corresponds to the slab geometry with diffusive walls, Fo from the first row of Table [III and 
other parameters from the left side of Table |l] Other lines differ from it as follows. The orange lines use Fo from the second 
row of Table [ll] The dashed red and orange lines use the right side of Table |l] The dashed blue line has absorbing container 
walls. The dashed pink line has the cylindrical container of radius b — 8a. The dashed black line is the hydrodynamic result 
in unlimited fluid. The vertical black lines give the effects 1-4, see text. 

exact treatment of the collision term, including higher Fermi-liquid parameters like F2, and allowing more general 
boundary conditions where degree of specularity depends on the angle of an incident quasiparticle. It seems unlikely, 
though, that these could lead to much better understanding of the experiments. 

IX. CONCLUSION 

We have calculated the response of Fermi-Bose liquid to an oscillating cylinder. This is applied to vibrating wires 
in •^He-^He mixtures. The results differ from ideal-gas results because of Fermi-liquid effects. In particular, the 
resonance frequency in the ballistic limit exceeds the ideal-fluid value because of "^He bound to quasiparticles and the 
Landau force, i.e. the elasticity of the Fermi liquid caused by interactions between the quasiparticles. The results 
are compared to measurements. For the comparison it is essential to take into account the size and the form of the 
experimental container. Good agreement is achieved without any fitting parameters. It seems that to explain the 
remaining differences would require 3D simulation of experimental volume. 
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FIG. 10: (a) The experimental (circles) and numerical results (solid lines) for wire 1 for two concentrations, 5.6% at saturated 
vapor pressure (green) and 7.0% at P = 10 atm (red). The hydrodynamic results (dashed lines) are indistinguishable on this 
scale. We see that the numerical results differ only slightly from each other in the ballistic limit. The reason for the difference 
between the experimental results is unknown, (b) The frequency parameters for wire 1, — 1.8% curve (green) and wire 2, 
xs — 9.5% curve (blue) are nearly equal, ~ 0.022. The numerical curves (solid lines) overlap in the hydrodynamic region, but 
start to differ around the bend in the curve, due to different values for Fi . The slopes of the calculated curves in the ballistic 
regime differ because of different values for Fq. The theoretical curves are identical to the principal curves in Fig. [9] 
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